Magnification Bias Corrections to Galaxy-Lensing Cross-Correlations 
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Galaxy-galaxy or galaxy-quasar lensing can provide important information on the mass distribu- 
tion in the universe. It consists of correlating the lensing signal (either shear or magnification) of 
a background galaxy/quasar sample with the number density of a foreground galaxy sample. How- 
ever, the foreground galaxy density is inevitably altered by the magnification bias due to the mass 
between the foreground and the observer, leading to a correction to the observed galaxy-lensing 
signal. The aim of this paper is to quantify this correction. The single most important determining 
factor is the foreground redshift zf. the correction is small if the foreground galaxies are at low 
redshifts but can become non-negligible for sufficiently high redshifts. For instance, we find that 
for the multipole I = 1000, the correction is above 1% x (5s/ — 2)/6/ for Zf > 0.37, and above 
5% x (5s f — 2)/bf for Zf > 0.67, where s/ is the number count slope of the foreground sample, 
and bf its galaxy bias. These considerations are particularly important for geometrical measures, 
such as the Jain and Taylor ratio or its generalization by Zhang et al. Assuming (5s/ — 2)/fe/ = 1, 
we find that the foreground redshift should be limited to Z/ < 0.45 in order to avoid biasing the 
inferred dark energy equation of state w by more than 5%, and that even for a low foreground 
redshift (< 0.45), the background samples must be well separated from the foreground to avoid 
incurring a bias of similar magnitude. Lastly, we briefly comment on the possibility of obtaining 
these geometrical measures without using galaxy shapes, using instead magnification bias itself. 

PACS numbers: 98.80.-k; 98.80.Es; 98.65.Dx; 95.35.+d 



I. INTRODUCTION 



The gravitational lensing of distant objects such as 
galaxies or quasars (39j located around a redshift Zb con- 
sists of both shear (a change in shape) and magnification 
(a change in size and flux), and is due to the integrated 
mass along the line of sight (LOS) between these back- 
ground galaxies and the observer. Correlating this lens- 
ing signal (cither shear or magnification) with the fore- 
ground galaxy number density at Zf < Zf,, we pick up the 
only part of the lensing mass which is correlated with the 
galaxy density at Zf, i.e. the matter located around Zf. 
Galaxy-lensing correlations arc therefore very interesting 
tools to study the matter distribution around the lenses. 
These correlations have been observed experimentally 
and studied theoretically both in the case of galaxy-shear 
correlation (also often called galaxy-galaxy correlations) 
1 1 2l jjjj- M. a, Ld, lal and galaxy-magnification correlation 

a, ESniil, liaLiS C3, LH, El, E3, El, til, [13, HH, H3, m, [ip . 

In addition to providing information on the matter dis- 
tribution, galaxy-lensing cross-correlations can be com- 
bined in order to build quantities which depend only on 
geometrical distances to the background and foreground 
samples, allowing measurements of dark energy (DE) pa- 
rameters independently of the exact expression of the 
matter power spectrum. These methods, that we call ge- 
ometrical methods, have been presented in (25l . [26l . [27j in 
the case of galaxy-shear correlations and in a more gen- 
eral way in [28[ for both galaxy-shear and shear-shear 
correlations. 



In this article, we would like to draw attention to 
the fact that these galaxy-lensing correlations depend 
crucially on estimates of the foreground galaxy density, 
which is inevitably altered by magnification bias, leading 
potentially to a systematic source of error. Magnification 
bias is a well known effect (a derivation can be found in 
H^l). Corrections due to magnification bias caused by 
the presence of the mass along the LOS between Zf and 
z = enter as follows. Suppose there is a mass over- 
density along the LOS. On the one hand, the observed 
number of galaxies per unit area on the sky is smaller due 
to the stretching of the apparent inter-galaxy spacing; on 
the other hand, very faint galaxies which were under the 
threshold of detectability are now visible, leading to a 
higher number density. The net effect depends on the 
number count slope s f and lead to a correction term be- 
ing added to the intrinsic galaxy density 



Sn = 5q+ Sn 



(1) 



where 5^ oc (5s j — 2). As a consequence, the galaxy- 
lensing cross-correlation gets an extra term due to the 
correlation between the lensing signal of the background 
galaxies and the extra foreground galaxy density 5^. How 
important this effect is, and under what condition this 
should be taken into account, is addressed as follows. 

In section UH we present the expressions for the mag- 
nification bias and for the galaxy-lensing correlations, in- 
cluding both galaxy-shear and galaxy-magnification cor- 
relations • Then, we study in section IIIII how the 
correction due to magnification bias depends on various 
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parameters in the problem: the background redshift Zb, 
the foreground redshift Zf, the angular scale, the selec- 
tion function width er, the number count slope Sf and 
the galaxy bias bf. In section ITVl we point out that mag- 
nification bias corrections are particularly important for 
the geometrical methods of [25| an d [23|- We discuss 
strategies for minimizing them. We also suggest a sim- 
ple extension of these geometrical methods which uses 
magnification bias itself, without the need for any galaxy 
shape information. We conclude in section[V] One recur- 
ring theme that would come up in our discussion is that 
the foreground redshift Zf is an important determining 
factor in the size of the magnification bias corrections. 
This is not only because a larger Zf makes for a larger 
lensing efficiency, it is also because a larger Zf gener- 
ally means intrinsically brighter galaxies and therefore 
steeper number count slope Sf. The latter is particulary 
true if precise (i.e. spectroscopic) redshifts are necessary 
for the foreground galaxies, such as in some versions of 
the geometrical methods. 



II. FORMALISM AND CORRELATION 
FUNCTIONS 

A. Magnification bias 

When magnification bias is taken into account, the ob- 
served galaxy overdensity at a mean redshift zq in the 
direction 9 is 



6„ (6, Zq) = S g (0, Zq) + 6^(8, Zq) 



(2) 



The first term is the intrinsic galaxy fluctuation inte- 
grated over a normalized window function W(z, z ) and 
is given by 



5 g {8,z ) 



dzW(z,z )b(z)d[x(z)d,z] (3) 



where x( z ) is the comoving distance to redshift z and 
S = (/> — p)j p is the matter overdensity. We consider 
here the case of a flat universe, but the expression can be 
easily generalized to an open or closed universe. We as- 
sume a scale independent bias such that S g = b(z)5; this 
should be correct for k < 0.05/iMpc -1 , though correc- 
tions could become important at smaller scale (see e.g. 
[30|). If the bias b varies slowly across the width of the 
selection function W, we can further assume that the 
bias is constant through the galaxy sample, allowing the 
simplification 

/>OC 

S g (6, z ) = b(z ) / dz W(z, z )6[ X (z)8, z] ; (4) 
Jo 

we will assume the relation holds in this paper, although 
it is straightforward to generalize our results to encom- 
pass both scale dependent and rapidly evolving bias. 



The second term of Eq. ([2]) encodes the overdensity due 
to magnification bias, and reads (assuming that s(z) is 
slowly varying across W(z,Zq)): 



6 M (0,zo) = (5a{zo)-2) I dz 

/o H{z) 



g(z,z )V 2 j_<f>[x(z)8,z] 
(5) 

where V\(f> is the 2D Laplacian of the gravitational po- 
tential in the plane perpendicular to the line of sight and 
s is the slope of the number count function (see e.g. [29j). 
For a survey with limiting magnitude m this is 



rflog 10 N{< m) 
dm 



(6) 



We have also introduced the lensing weight function 

jX{z')-x{z), 



9(z,z ) = x{z) I dz 



X(z') 



-W(z',z ) (7) 



that can be thought of as proportional to the probability 
for a source around z to be lensed by mass at redshift z; 
it is peaked at a redshift corresponding roughly to half 
the comoving distance to the source (for a review of weak 
lensing, see e.g. (3lj). 

If we make use of the Poisson's equation, Eq.© be- 
comes 



6ti(0, zo) 



3H? 
2 c 2 



I Q m ,o) (5s(z ) - 2) 



(8) 



x / dz — — 

lo H(z) 



g(z,z )(l + z)6\x(z)8,z]. 



B. Cross-correlations 



Galaxy-shear correlations 



Galaxy-shear correlation (also often called galaxy- 
galaxy lensing) has been the object of considerable in- 
terest from theoretical and observational points of view 
0, 0, H, H, H, @, 0, 11] ■ It consists of correlating the cosmic 
shear of background galaxies located around a redshift z b 
with the number density of foreground galaxies at red- 
shift Zf. Here, we focus on the electric-component of the 
shear (sometimes called the scalar component or the gra- 
dient mode) and we will use n to refer to this scalar part: 
the convergence. The cosmic shear of the sources is due 
to all the mass between z b and z = 0; however correlating 
this shear with galaxy number density at Zf, we extract 
the part of the lensing mass that is correlated with the 
galaxy density at Zf, i.e. the matter located around Zf. 
Therefore, galaxy-shear correlations allow us to get im- 
portant information about the galaxy-mass power spec- 
trum. 

More precisely, due to magnification bias, the galaxy- 
shear two-point correlation function is the sum of two 
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terms 

w nK (0,z f ,z b ) = (6 n (6,Zf)K(d',z b )) 

= (5 g (0, z f ) K (0', z h )) + (5 M (0, Zf)n{tf, z b )) 
= w gK (9, Zf, Zb) + Wp K (0, Zf, Zb) 

(9) 

where cos 9 = 0.0'. We then consider the Legendre co- 
efficients C 9 t K {zf, Zb), Cg K (zf, Zb) of the correlation func- 
tions. They are defined as 



gK (9,z f ,z b ) = 



2£+l 

47T 



Cf K (z f ,z b )P e {cosO) (10) 



where the Pi are the Legendre polynomials. The galaxy- 
shear angular power spectrum will then also be the sum 
of two terms, 

cr(z f , z b ) = ci K ( Zf ,zb) + cr(z f , z b ). (ii) 

We calculate these two angular power spectra Cf K , C^ K 
using the Limber approximation which is accurate for 
large multipoles I > 10 [32]: 



Cf K (z f ,z b ) 



1 3 

2 c 



2 Q m ,o ) bf 



° d W(z,z f ) 9 (z,z b ) 
X (z) 



3H, 



Cjr(zf,z b )=[--§-Sl m ,o) (5s, -2) 



(12) 



2 c? 

(13) 

where P(k, z) is the 3D matter power spectrum for a wave 
number k and at a redshift z, bf the intrinsic clustering 
bias of the foreground galaxy and s / the mean number 
count slope of the foreground galaxy sample. Usually, 
the signal we are interested in is C| K , while C^ K will 
constitute the magnification bias induced correction to 
this signal. It is then important to know to what extent 
this correction can be neglected or should be taken into 
account. We will discuss these issues in detail in the 
following sections, but before doing so, we would like 
to emphasize that the same sort of corrections exist for 
galaxy-magnification cross-correlations. 



2. Galaxy-magnification correlations 



ic , im, yj, 

24l | is very 



Galaxy-magnification cross-correlation [i 

M M \m M Ei M M M, M 11, m 

similar to galaxy-shear correlation replacing the shear by 
the magnification of the background galaxies. In order 
to measure these correlations, one should correlate the 
galaxy densitiy of a background distribution with the 



density of a foreground distribution leading to the corre- 
lation function 



(14) 



w nn {9,z f ,z b ) =w gg (9,Zf,z b ) +w gfJl (9,Zf,z b ) 

+ w f _ ig (9, z f , z b ) + 11)^(6, Zf, z b ) 

If the background and foreground distributions do not 
overlap, and they are far apart, it is easy to see that the 
intrinsic galaxy-galaxy correlation w gg (9, Zf, z b ) and the 
correlation between magnification of foreground galaxies 
and background galaxy density w^ g (9, Zf, z b ) both van- 
ish. The two remaining terms 

w nn {9,zf,z b ) = w gf j,(6,Zf,z b ) + w^(9,zf,z b ) (15) 

are comparable to Eq.@. 

Similarly to C™ K , the angular power spectrum will be 
the sum of two terms 



C^(z f ,z b ) = C 9 /(zf,z b ) + C^(z f ,z b ), (16) 



where 



Ct(z f ,z b ) 



3 #o 2 



^m,o bf(5s b - 2) 



2 c 2 

d W(z,Zf) 9 (z,z b ) 
o X 2 (z) 



(17) 



(18) 

We can see that these expressions depend strongly on 
the galaxy samples, through the bias bf and the number 
count slopes Sf,s b . Here, Cjf constitutes the magnifica- 
tion bias correction to C| M , the quantity that observers 
hope to measure by cross-correlating number counts of 
two widely separated samples of galaxies. 



3. Formalism that combines the two types of correlations 



If we compare Eq. (|12H13j) and Eq. (|17H18j) . it is easy to 
see that they have a similar structure. To give a uni- 
fying treatment, let's introduce the matter-convergence 
and convergence-convergence power spectra 

Cf K (zf,z b ) = (^-^-^m,o 

x r d z W ^%]f Zb \ l + z)P(l/ X (z),z), 
Jo X (Z) 



(19) 



G i K (zf,Z h ) = (^-jf^m-O 

Jo H (z) x 2 {z) 



(20) 
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then, we can relate the different power spectra in the 
following way 



Cf K (zf,z b ) 
C>f K (z f ,z b ) 



b f Cf K (z f ,z b ) (21) 
(5s f -2)Cr(z f ,z b ) (22) 



and 



Ct(z f ,z b ) 



(5s b -2)b f Cr(z f ,z b ) 
(5s b -2)(5sf-2)CH K (z f ,z b 



(23) 
(24) 



We are interested in the fractional magnification correc- 
tions to the observed galaxy-lcnsing correlations, and 
in both the galaxy-shear correlation and the galaxy- 
magnification correlation, that ratio is 



C$ K {zj,z b ) _ C^(z f ,z b ) _ (5a/ -2) 



Cr(z f ,z b ) Cr(z f ,z b ) 



where r is defined as 



C£ K (z f ,z b ) 



r(z f ,z b ;£) (25) 



(26) 



In what follows, we will display several figures showing 
r, but it should be kept in mind that ultimately, the 
fractional correction of interest is r multiplied by (5s f — 
2)/bf, a quantity which is sample dependent. It suffices 
to say that for vast classes of galaxy samples, this factor 
of (5s f — 2)/bf could be easily of order unity. We will 
discuss this point in detail in i jlHDI 

In all illustrative examples below, we will consider for 
simplicity a gaussian selection function of width a 



W(z,z ) 



1 



2ira 



cxp 



2ct 2 



(27) 



and the following cosmological parameters: the Hubble 
constant h = 0.7, matter density fi TO = 0.27, cosmo- 
logical constant Qa = 0.73, baryon density Q b = 0.046, 
power spectrum slope n = 0.95 and normalization erg = 
0.8. We employ the transfer function of [33[ and the pre- 
scription of |34| for the nonlinear power spectrum. 



III. CORRECTION DUE TO MAGNIFICATION 
BIAS 

We would like to understand better to what extent 
the correction due to magnification bias should be taken 
into account, and if so, in what cases. Before going into 
details, we would like to emphasize that the correction 
to Cf K due to C<t K can be non-negligible; in Fig. [T] we 
plotted the correlation "signal" Cf K (zf, z b ) along with 
the corrected observable Cf K (zf, z b ) + C£ K (zf, z b ) for two 
foreground and background redshifts. We see that the 
correction is very small for low foreground redshifts, but 
this is not the case for a higher Zf. 

In the following, we will study in detail the ratio in 
Eq. (|2l)|) . or its cousin r in Eq. (|2"6")) . and their dependence 
on the different parameters. 




1000 



FIG. 1: Galaxy-lensing power spectrum Cf K (zf,zt) and the 
corrected power spectrum Cf K (zf,zt) + C( K (zf, z b ) for two 
redshift configurations. The redshift width of the selection 
function is a = 0.07. 



A. Dependence on the background galaxies 
redshift 25 

The ratio r depends on z b through the lensing weight 
function g. In Fig. ® we can see that r grows when the 
background sample gets closer to the foreground galax- 
ies. This is easy to understand: the galaxy-lcnsing corre- 
lation, which is proportional to C\ K , becomes small when 
the source-lens separation becomes small. More quanti- 
tatively, if one approximates the selection function by a 
delta function i.e. W(z, zf) ks S(z — zf (which is a rea- 
sonnable approximation if the width a of the selection 
functions is small and if the sources and lenses distribu- 
tions do not overlap), we have: 



Ct(z f ,z b ) 



— — 5-i 'm,0 



X(l 



X(zb) ~ x(zf 
X(zfx(zb) 
z f )P(£/x(zf,zf 



(28) 



showing that Cf K goes to zero for close background and 
foreground samples; this is due to the fact that when the 
lenses are very close to the sources, the lensing weight 
function g is almost null (it would exactly vanish if in- 
deed W were a perfect delta function). Under the same 
approximation, we have 



cr(z f ,z b )^(J^-n m ,o 



dz — ^— 

H(z) 



X(zf 



X(zb) 



(29) 



5 




0.001 



0.0008 



FIG. 2: r(z f ,z b ;£) = C? K (z f , z b )/Cf K {z f , z b ) for two differ- 
ent fixed foreground redshifts Zf = 0.1, 1 and for three values 
of £ = 10, 100, 1000. The width of the selection function is 
a — 0.07. For each zj we let z b vary in [Zf,Zf + 1]. We can 
see that the ratio r decreases with the separation between the 
background and foreground samples. 



will explore this further next. 



B. Dependence on the foreground galaxies redshift 




which tends toward a finite value when Zb — * Zf. It 
is therefore not surprising to see that the ratio r be- 
comes large when the background galaxy redshift ap- 
proaches that of the foreground. Of course, for real- 
istic observations, the selection function is never quite 
a delta function, and the precise behavior at low back- 
ground redshift will depend on its precise shape. We will 
study that aspect in more detail in ^III CI (see also HH ) . 
Note also that when Zh approaches z/, the whole premise 
of galaxy-magnification cross-correlation measurements 
breaks down: as described in i jHB 21 the idea of such 
measurements is to correlate the number counts of a fore- 
ground sample and a background sample which are suffi- 
ciently far apart so that the intrinsic clustering term w gg 
is unimportant. 

Another interesting feature of Fig.[2]is that for high zn, 
r asymptotes to some minimal value (independent of the 
background position) . This can also be understood easily 
through the equations (f2"5|) . (pZTJ)) ; indeed, when x( z f) ^ 
x{zb), x( z b) cancels out of Eq. ([28|) while Eq. (j29|) becomes 
weakly dependent on x( z b)- Therefore, if one's goal is to 
minimize (magnification bias induced) corrections to the 
galaxy- lensing correlations, it is important to make sure 
the sources and the lenses are widely separated. However, 
it is worth noting that according to Fig. O r can still be 
fairly substantial even when Zh differs significantly from 
zj. This holds especially if zt is sufficiently large. We 



FIG. 3: Ratio r(zf,z b ;£) for different foreground redshifts zf, 
for each 2/, we choose a background redshift z b = Zf + 1 , such 
that according to Fig. [2] r is minimized. The redshift width 
of the selection function is a = 0.07. On the top panel, low 
multipoles are shown, while the bottom plot shows r for high 
multipoles. 



The foreground redshift appears to be the crucial pa- 
rameter which determines when the correction to C| K due 
to Cf K is large. Indeed, the correlation Cf K comes from 
the lensing of both background and foreground galaxies 
by the mass located between Zf and z = 0. The lens- 
ing efficiency associated with this mass increases as zj 
increases. This effect is already visible in Fig. [2] but is 
even more striking in Fig. [3] in which we plot r as a func- 
tion of Zf (for each point, we choose Zb = Zf + 1 so that 
r is close to its asymptotic minimum value for this Zf). 
Obviously, the curves for low multipoles have a very dif- 
ferent behavior from the curves for high multipoles (we 
will study in more detail this dependence on £ below). 
However, the main point is that for each multipole £, r 
grows with z f and reaches a non-negligible value for high 
Zf. For example, for £ = 1000, r = C£ K /Cf K is above 1% 
as soon as Zf > 0.37, and above 5% for Zf > 0.67. As 
a consequence, magnification bias correction should be 
taken into account when interpreting galaxy-lensing cor- 
relation measurements if the foreground redshift is suffi- 
ciently high. 

This dependence on the foreground redshift Zf also al- 



6 



lows us to infer that this correction is negligible for the 
cross-correlations that have been observed until now. In- 
deed, in the case of the recently observed galaxy-shear 
correlations (or so called galaxy-galaxy lensing) using 
the SDSS survey S 1J3, I|, the lenses (fore- 
ground) are mostly distributed [l| around a mean rcdshift 
(zf) = 0.1. It is safe to assume that these galaxies have 
|5s/ — 2\/bf of at best unity Therefore the fractional 
correction is small: \C^ / Cf h \ < 1%. Similar conclu- 
sions hold in the case of galaxy-magnification correlations 
measured by cross-correlating foreground galaxy counts 
with background quasar counts, as carried out by [2l| : 
here again, the distribution of the foreground galaxies is 
peaked around the low redshift (zf) = 0.24, leading to 
a less than 5% correction even for a large number count 
slope. 

The situation will be quite different for future weak 
lensing surveys, like LSST [H or DUNE (now EUCLID) 
[37l ]; indeed, these surveys will obtain precise photomet- 
ric redshifts up to very high redshifts z ~ 3, and therefore 
the magnification bias correction will not be a priori neg- 
ligible. In particular, for high lens redshift Zf > 1, if the 
foreground number slope Sf is not such that the magni- 
fication is small, the correction can be of order one. To 
keep the magnification bias correction small, one can ei- 
ther focus on low Zf galaxies, or select samples where s/ 
is close to 0.4 (or weigh galaxies by l/(5s/ — 2)). Perhaps 
a better strategy is to include the magnification bias cor- 
rection from the beginning when interpreting measure- 
ments. 



2 in the case of non-linear matter power spectum (case 1); 
in order to emphasize the effect of non-linearities, we also 
plot r(i) assuming the matter power spectrum is equal 
to the linear one (case 2). (All figures in the rest of this 
paper use the nonlinear power spectrum.) 

We can see that r generally decreases with I. We can 
also see a bump around I ~ 1000 in the non-linear case. 
Since non-linearities become important at different £'s for 
Cg K and Cf K (because the former is sensitive to fluctu- 
ations between redshift zero and Zf while the latter is 
sensitive to fluctuations at around Zf), it is not surpris- 
ing that their ratio is not monotonia At sufficiently high 
^'s, the ratio r asymptotes to some low value. 

The overall behavior of r can be understood as fol- 
lows. Approximating once again the selection function 
as a delta function, we need only to consider the ratio of 
Eq.([29|) to Eq.([28]). At sufficiently high £'s, the integral 
over z in Eq. (|29|) is such that it is dominated by z's where 
P(£/x(z), z) is more or less power-law in I i.e. P oc £~ 3 
(originating from the fact that P{k) ex k~ 3 at high fc's). 
Eq. ([28|) is likewise roughly proportional to £~ 3 for suffi- 
ciently high C&. Their ratio r is therefore independent 
of £, giving rise to the high £ asymptote that we see in 
Fig. |4] At low £'s, on the other hand, the integral in 
Eq. ([29]) goes through the maximum of the matter power 
spectrum, and the result depends strongly on £, leading 
to the steep decreasing slope for low £ noticed in Fig. 0] 

We should emphasize that at high £'s (> I0 3 ), the 
actual ratios of interest Cf K /Cf K and Cf^/C^ (Eqs. 
[12] . [13"] . [17] . [TB] ) are affected by nonlinear and scale 
dependent galaxy bias. 



Dependence on the multipole £ and the 
selection function width a 



case 1 : NL Power Spectrum 
case 2: Linear Power Spectrum 




10' 



10" 



Multipole 1 



10" 



10 



FIG. 4: Ratio r(z f , z b ;£) = C^{z f ,z b )/Cl K {z s , Zb) as a func- 
tion of the multipole I for Zf = l,Zt = 2 in two cases. Case 
1: r is computed using a non-linear matter power spectum 
[34| ; Case 2: r is computed assuming the matter power spec- 
trum is equal to the linear one. For low £, r decreases steeply, 
before reaching a roughly constant plateau for high £. 



As discussed above, the correction due to C" K depends 
on £\ in Fig. [5] we plot rasa function of £ for Zf = 1, Zf, = 



1.4 
1.2 \ 
1 

0.8 

0.6 

0.4 

0.2 


1 



a=0.07, 1=1000 
a=0.15, 1=1000 
a=0.30, 1=1000 



1.2 



1.4 1.6 



1.8 



FIG. 5: r(zf, zi,; £) = Cf K (zf, zt,)/C e K '(zf , Zb) as a function of 
Zf, for Zf = 1 and I = 1000; three different widths of the selec- 
tion function are shown: a = 0.07, 0.15, 0.30. For background 
redshifts close to foreground ones, r depends strongly on a, 
decreasing as a increases; for high z b > z/ + 3cr, the width of 
selection function has no influence anymore, as expected from 
the delta function approximation. 



In this paper, we are mostly interested in the cases 
where the foreground and background distributions do 
not overlap, and for which the width a of the selection 
function is of little importance, as long as it is small (the 
case where the foreground and background galaxies are at 
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the same redshift has been studied extensively in [35j] ) . 
This is confirmed by Fig. [5] where we see that as long 
as the separation between foreground and background 
is large compared to er, the precise value of a does not 
matter much. In cases where z b is actually close to zf, the 
main effect of increasing a is to make Cf K larger, hence 
making the ratio r smaller (recall that Cf K vanishes if 
the selection function is an exact delta function). 



simple way to implement these techniques using galaxy 
counts alone, without using any galaxy shape informa- 
tion. Second, magnification bias corrections may affect 
these measures in a significant way, or more precisely, af- 
fect significantly the inferrence of dark energy parameters 
from these measures. 



D. Dependence on the bias bf and number count 
slope Sf 

In the previous sections , we have studied the pa- 
rameter dependence of the ratio r = C£ K /Cf K . How- 
ever, the quantity of observational interest is actually 
^ 5s ^ ^ r (Eq.(|25p). which depends crucially on the lens 

sample dependent quantities (5s / — 2) and bf. The ratio 
(5s f — 2)/bf can vary anywhere from —1 to 2 depend- 
ing on the sample limiting magnitude and redshift, and 
it generally increases as one focuses on brighter galax- 
ies and higher redshifts (see Fig. 1 of [1^]). It is worth 
noting that measurements that require spectroscopic red- 
shifts (such as in the geometrical method proposed by 
[25l ] . discussed below) focus on bright galaxies for which 
(5s/ — 2)/bf is generally large i.e. order unity or above. 



IV. GEOMETRICAL METHODS AND 
MAGNIFICATION BIAS 



0.705 



^ 0.704 




FIG. 6: Dependence of the ratio R = 
on the DE equation of state w for zj 
Note that R is a weak function of w 



Cf 



(z f , zi )/cr(z f ,z 2 ) 

0.3, zi = 0.5, z 2 = 0.7. 



A. Extensions using number counts 

Jain and Taylor [25| (see also [H, [2?]]) proposed using 
the ratio of galaxy-shear correlations (i.e. in a galaxy- 
galaxy lensing experiment, as described in i )IIB ip to 
constrain a certain combination of angular diameter dis- 
tances: 



R = 



Cf K (z f ,z 2 ) (x(z 2 )-x(zf))x(zi) 



(30) 



The second equality holds if the selection function are 
sharply peaked; this places an especially high demand on 
the redshift accuracy of the foreground galaxies (2?| ■ 

It is fairly straightforward to see that the same ratio 
can be obtained by galaxy-magnification correlation mea- 
surements, the type described in miB 21 In other words, 
by cross-correlating the number counts of a foreground 
sample and a background sample of galaxies, one can iso- 
late Cf^ (Eq. |I7| , temporarily ignoring the magnification 
bias correction that is in Eq. |I8| ) . If one can measure this 
cross-correlation using two different background samples 
but with the same foreground, we can form a ratio anal- 
gous to the Jain and Taylor ratio of Eq. (|30l) : 



Cr(z fl zj) = (5si - 2) ( X (zi) - X (zf)) X{Z2) 
' Cr(z f ,z 2 ) (5s 2 -2)( X (z 2 )- X (z / )) X (z 1 ) 

(31) 

where the second equality holds again under the assump- 
tion of a sufficiently narrow selection function. This ra- 
tio differs from Eq. (|3"0"|) only by factors that involves the 
number count slope which can be independently mea- 
sured. 

Similar logic applies to generalizations of the Jain and 
Taylor ratio, such as the one proposed by Zhang et al. 
[281 ] . It exploits what is known as an offset linear scaling 
of the galaxy-shear correlation (similar scaling applies to 
shear-shear correlation as well): 



By "geometrical methods", we refer to techniques that 
involve combining lensing measurements at different red- 
shifts to obtain measures that are purely geometrical in 
nature (i.e. they depend only on combinations of geo- 
metrical distances, and not on the power spectrum or its 
growth). There are two simple observations we would 
like to make about these techniques. First, there is a 



Cf K (z f , z b ) = F e {z f ) + G e (z f )/ Xc s(z b ) (32) 



where 



1 



Xes(z b ) 



dz 



, W(z',z b ) 
X(z') 



(33) 
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and 



Fi(z f ) = (Iff-O^b/ 



Gt{z f ) 



X ' dz W{z : Z f\ l + z)P(t/ X {z),z) (34) 
'o X\ z ) 

dzW(z,z f )(l + z)P(£/x(z),z). (35) 



Zhang et al. [28| suggested combinations of C^ K 's from 
different samples which are purely geometrical in nature, 
such as 

R z = Cf K {z f ,z 3 ) - Cf K {z f ,zi) = Xcff(z 3 )~ 1 -Xcff(zi)~ 1 

" C r{ z fi z 2)-Cf K {zf,Zi) Xcff(z2) _1 -Xoff(2l)" 1 

_ (36) 

An advantage of this approach compared to [251 ] is that 
it is less demanding on the redshift accuracy. 

Once again, one can implement the same idea using 
the cross-correlation of number counts alone i.e. one can 
study instead: 

_ Cf (z f , z 3 )/(5s 3 - 2) - Cf {z f ,zx)l{hs x - 2) 



C 9 ^(z f , z 2 )/(5s 2 - 2) - Cf" (z/,zi)/(5 Sl - 2) 

Xeff^r 1 - Xcff(^l) _1 



Xes(z 2 ) 1 - Xoff(^l)" 



B. Magnification bias corrections 



(37) 




FIG. 7: Bias in the inferred dark energy equation of state 
Aw if magnification bias corrections are ignored in inter- 
preting the Jain and Taylor ratio. This is as a function 
of the background redshift z\, for three foreground redshifts 
Zf = 0.15, 0.35, 0.45. For each z\, we fix the other background 
redshift zi — z\ + 0.4 (our results do not depend sensitively 
on the choice of Z2 — z\). The selection functions are approx- 
imated by delta functions, and (5s/ — 2)/bf ~ 1 is assumed. 
The calculation is done for £ = 1000. 



The methods described above are interesting because 
they allow us to isolate from cross-correlation measure- 
ments purely geometrical constraints on dark energy. 



Such constraints can be compared with constraints from 
the growth of structure as a consistency test of general 
relativity. Fig. [5] shows the Jain and Taylor ratio R as 
a function of the dark energy equation of state w. We 
would like to draw attention to the fact that R is a weak 
function of w. This means an accurate determination of 
w would require a very precise measurement of R. For 
instance, in the configuration of Fig. [5J to constrain w to 
5%, we would need to measure R to within ~ 0.04%. This 
suggests the magnification bias corrections we have been 
studying might be especially relevant: their presence, if 
unaccounted for, would lead to incorrect inferrence on 
the dark energy equation of state. 

More quantitatively, the observed Jain and Taylor ratio 



Robs 



cr(zf,z 2 )+cr(z f ,z 2 ) 

C!"(z f , Zl ) ( l+ { -^r(z f ,z i; e) ' 
Ct K (zf,z 2 ) 1 r{zfjZ2 . £) 



(38) 



where 2/ is the foreground redshift, and z\ and z 2 are two 
different background redshifts. Inferring w from i? bs, 
without accounting for its difference from the idealized 
ratio R in Eq. (|38|) (i.e. using R as the model for the 
data), leads then to an inferred DE equation of state 
^inferred which differs from the true one ui trU c; we can 
Aw = u>true — ^inferred the resulting bias. In Fig. [71 
this bias Aw is shown for various redshift configurations 
and I = 1000. Here, we consider redshifts ranging up 
to z = 2.6, approximate the selection functions by delta 
functions, and assume (5s/ — 2)/6/ ~ 1. The approxima- 
tion and assumption are reasonable given that the Jain 
and Taylor method requires essentially spectroscopic red- 
shifts for the foreground galaxies: such galaxies are gen- 
erally bright and their values for (5s/ — 2)/6/ are typically 
large (see [29[). We find that for very low foreground red- 
shifts (zf < 0.15), the bias Aw is generally less than 5% 
for all combinations of the two background redshifts z\ 
and z 2 . For intermediate redshifts 0.15 < Zf < 0.45, 
whether one can keep Aw below 5% depends on the 
Zi,z 2 configuration: for low z±,z 2 (i.e. approaching Zf), 
the magnification bias correction is simply too large to 
be ignored; for sufficiently high Zx,z 2 , r(zf,zi;£) and 
r(zf,z 2 ;£) becomes close to each other, almost indepen- 
dent of the background redshifts (see Fig. [5]) and so the 
magnification bias corrections cancel each other to a large 
extent in i? b s - Finally, for high foreground redshifts 
(zf > 0.45), the magnification bias corrections are simply 
too large to allow such delicate cancelations: Aw exceeds 
5% for all background configurations. 

One could of course, when interpreting data, account 
for the magnification bias corrections from the beginning 
and avoid such a bias Aw. Note however that these cor- 
rections depend on the power spectrum, making the ob- 
served ratio i?obs not strictly geometrical anymore. In 
other words, this ratio is a useful, pure geometrical mca- 
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sure only in so far as the magnification bias corrections 
are small. 



V. CONCLUSION 

In this paper, we present a detailed study of the correc- 
tion due to the magnification bias of foreground galaxies 
in the context of galaxy- lensing cross-correlations. We 
find that the most important factor that determines the 
size of this correction is the foreground rcdshift: the size 
increases with rcdshift. For example, we found that for 
t = 1000, the correction is above 1% x (5s f — 2)/bf for 
z f > 0.37, and above 5% x (5s f - 2)/b f for z f > 0.67. 
Maintaining a low foreground rcdshift is therefore neces- 
sary to keep the magnification bias correction small. (We 
find a sufnccntly large separation between foreground and 
background, Az ~ 1, is also necessary.) Or else, one 
should explicitly account for the magnification bias cor- 
rection when interpreting data. This is especially true for 
future surveys which are expected to reach high precision 
and high redshifts. Finally, we present a simple extension 
of existing geometrical methods to measure dark energy 
parameters independently of the matter power spectrum. 
In the case of galaxy-shear [25|, [28| and shear-shear (28| 



correlation functions: the galaxy-shear correlation can 
be replaced by galaxy-magnification correlation, allow- 
ing measurements that involve only galaxy counts, with 
no galaxy shape information necessary. However, we 
also point out that these geometrical methods are gener- 
ically sensitive to magnification bias corrections. For 
(5s f — 2)/bf = 1, we find that the foreground rcdshift 
should be limited to Zf < 0.45 in order to avoid biasing 
the inferred dark energy equation of state w by more than 
5%, and that even for a low foreground redshift (< 0.45), 
the background samples must be well separated from the 
foreground to avoid incurring a bias of similar magnitude. 

As this paper was nearing completion, we became 
aware of a preprint by Bernstein [H| that has some over- 
lap with the magnification bias issues that we discuss 
here. 
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